%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Archivo parteE.m
% Calcula el tiempo medio por cliente en el sistema


rand('seed', 3230);
max_time = 300; % 5 horas de atencion

i=1;
[tiempos(1) xx] = simulate(1, max_time);
medias(1) = tiempos(1);

i = 2;
while (i <= 10)
	[tiempos(i) xx] = simulate(1, max_time);
	medias(i) = (medias(i-1) * (i-1) + tiempos(i)) / i;
	i = i + 1;
end

fprintf(stdout, 'Tiempo medio de cliente en el sistema para 10 simulaciones: %g\n',
		medias(i-1));
alpha = .05;
inv = norminv(1-alpha/2);
semi = inv*sqrt(var(medias(1:i-1))/(i-1));
fprintf(stdout, 'Semiamplitud del intervalo de confianza para alpha = %g: %g\n',
		alpha, semi);

title  'Tiempo medio de cliente en el sistema'
xlabel 'Iteracion'
ylabel 'Longitud de Cola de Recepción'

plot(1:10,tiempos);

tiempos'
